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ABSTRACT 

We obtained long-slit spectra of high signal-to-noise ratio of the galaxy M32 with the 
Gemini Mult i- Object Spectrograph at the Gemini-North telescope. We analysed the 
integrated spectra by means of full spectral fitting in order to extract the mixture of 
stellar populations that best represents its composite nature. Three different galactic 
radii were analysed, from the nuclear region out to 2 arcmin from the centre. This 
allows us to compare, for the first time, the results of integrated light spectroscopy 
with those of resolved colour- magnitude diagrams from the literature. As our main 
result, we propose that an ancient and an intermediate- age population co-exist in 
M32, and that the balance between these two populations change between the nucleus 
and outside one effective radius (Ireff) in the sense that the contribution from the 
intermediate population is larger at the nuclear region. We retrieve a smaller signal of 
a young population at all radii whose origin is unclear and may be a contamination 
from horizontal branch stars, such as the ones identified by Brown et al. in the nuclear 
region. We compare our metallicity distribution function for a region 1 to 2 arcmin 
from the centre to the one obtained with photometric data by Grillmair et al. Both 
distributions are broad, but our spectroscopically derived distribution has a significant 
component with [Z/Zq] ^ — 1, which is not found by Grillmair et al. 

Key words: stars: atmospheres, stars: evolution, galaxies: abundances, galaxies: 
evolution, galaxies: stellar content 



1 INTRODUCTION 

M32 is a controversial galaxy. This low-mass satellite of 
M31 has long been considered the prototyp e of the compac t 
elliptical morphological classification (Ben der et al.l Il992h: 
low mass, high surface brightness galaxies, tidally truncated 
companions to massive galaxies. So far, few objects match 
this description (Z iegler & Bender 1998; Mieske et al. 200 5; 
IChilingarian et al.ll2QQ7l ). The origins of the peculiar struc- 
tural properties of M32 are still a matter of debate and the 
proposed models span a wide range of hypotheses: from a 
true elliptical galaxy at the lower extre me of th e mass se- 
quence to a threshed early-type spiral (|Nieto PrugnieL 



* Based on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under a cooperative agreement with the NSF 
on behalf of the Gemini partnership: the National Science Foun- 
dation (United States), the Particle Physics and Astronomy Re- 
search Council (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research Council 
(Australia), CNPq (Brazil), and CONICET (Argentina) (observ- 
ing run ID: GN-2004B-Q-74). 



19871: iBekki et al.1 l200ll: IChoi et al.1 l20Q2l : iGrahamI l2002l : 
Gallagher et al.ll2005l : lYoung et al.ll2008l ). 



Controversial or not, M32 is a galaxy we need to un- 
derstand because it plays an important role in resolving 
controversies about how to interpret the integrated spec- 
tra of galaxies: it is an early- type galaxy where we can get 
extremely high signal-to- noise (S/N) in the spectra, good 
coverage outside the nucleus to look for radial population 
gradients and resolve individual stars to a very meaningful 
magnitude level. If we hope to claim that we understand the 
evolution of distant early- type galaxies, we must first be able 
to demonstrate an understanding of M32 stellar population. 

The stellar population of M32 has been extensively 
studied by means of integrated spectroscopy and re- 
solved photometry. Most studies of the integrated spec- 
trum of M32 have found signatures of an intermediate^ 
a^e a nd close to sol a r metallicity population (e.^. Ros 
1985; "Davidge" "l990'; Freedman 1992; Jones k Worthe 
1995; Vazdeki s & Arimoto 1999^ del Burgo et al. 2001; 
ISchiavon et al ] \200i : IWorthevll2004l : iRose et al.ll2005h . The 
most commonly used technique of deriving mean- ages and 
mean-met allicities f r om in tegrated spectroscopy, as artic- 
ulated by IWorthevI (|l994l ). consists in comparing selected 



2 P. CoelhOj C. Mendes de Oliveira and R. Cid Fernandes 

Table 1. Parameter coverage of the fits. Adopted: Zq = 0.017. 



Model 


log(age[yr]) 






[z/z©] 




Wavelength 


range (A) 


BC03 


8.0 - 10.3 


-2.2, 


-1.6, 


-0.6, -0.3, 


+0.05, +0.5 


4700 - 


8750 


Chariot & Bruzual 


8.0 - 10.3 


-2.2, 


-1.6, 


-0.6, -0.3, 


+0.05, +0.5 


4700 - 


7400 


PEGASE-HR 


8.0 - 10.3 


-1.6, 


-0.6, 


-0.3, +0.1, 


+0.05 


4700 - 


6800 


Vazdekis et al. 


8.0 - 10.3 


-1.6, 


-1.2, 


-0.6, -0.3, 


+0.05, +0.25 


4700 - 


7400 



spectral indices to the predictions of simple stellar popu- 
lations (SSP) models. Quoting an SSP-age and an SSP- 
metallicity, however, d oes not ne cessarily im ply that a 
galaxy is an SSP. In fact jBica et al.l (|l990. ) and .Hardv et alj 
(| 19941 ) have shown that the spectral energy dist ribution of 
M32 cannot be adequately fit by SSP models, and lBica et al.l 
([l990) argued that oversimpUfications might explai n some of 
the c o nflicting results foun d in literature (see also [Oconnel. 
198Q| : ISchmidt et al.1 Il989l : iMagris Bruzuall Il993l : IRosi 



1994l : lTrager et al.ll200Ql ~ 



Resolved phot o metric studies include t hose of 

lAlonso-Garcia et al.1 (|2004l l and I Srown et aP (|2000l l 
who found signatures of a true ancient population in M32 
[by the identification of RR-Lyrae and an extended hori- 
zontal branch (HB), respectively], and [Oavidge & Jensen 
(|2007l l and references therein, who identified signatures of 
an intermediate age population. The most comprehensive 
study to date of the resolved stellar populat ion of M32 has 
been carried out by iGrillmair et al.l (|l996h , who resolved 
individual stars down to slightly below the level of the 
HB with the Hubble Space Telescope (HST) Wide-Field 
Planetary Camera 2 (WFPC2) in a region 1—2 arcmin 
from the centre of the galaxy. The most striking feature 
of this colour-magnitud e diagram (CMD; and also the one 
bv IWorthev et alTl2004l , at a larger radii) is its composite 
nature, implying a range in [Z/H] from roughly solar down 
to below —1.0 dex. 



Spreads in age and Z are as important as the mean val- 
ues in understanding the star formation history of a galaxy, 
and a study that takes into account the composite stellar 
population of M32 may lead to a better understanding of its 
origin and evolution. In order to address this question, we 
obtained spectra of M32 along its major axis out to a radius 
of 2 arcmin, with the Gemini-North Telescope and Gem- 
ini Multi-Ob ject Spectrograph (GMOS ). We adopt the code 
STARLIGHT ^Cid Fernandes et alll2005ll and state-of-the-art 
stellar population spectral models teruz ual k Charlotl[2QQ3l . 
hereafter BC03; Le Borgne et al. 2004, hereafter PEGASE- 
HR; Chariot & Bruzual in preparation, hereafter Chariot & 
Bruzual; Vazdekis et al. in preparation, hereafter Vazdekis et 
al.) to perform a pixel- by-pixel modelling of the integrated 
spectrum of M32. For the first time we extend the analy- 
sis of integrate d spectra of M32 to th e radii of the resolved 
CMD study bv lCrillmair et al.1 (|l996l ). 



2 OBSERVATIONAL DATA 

Long-slit spectra of M32 were obtained on 2004 October 
10, 11 and 16 with the GMOS on Gemini-North under the 
Program GN-2004B-Q-74. We used the R400 grating, a 0.75- 
arcsec slitwidth and an effective target wavelength at 6800 A, 
resulting in spectra with coverage 4700 — 8930 A and an 
average resolution of full width at half-maximum (FWHM) 
5.4A. 

Four exposures of 30 s each were taken with the slit po- 
sitioned over the nucleus of the galaxy, along the major axis 
of the galaxy [position angle (PA) = 165°]. The slit was sub- 
sequently moved along the major axis off the ce ntral region 
to av oid contamination by scattered light (see [r ose et al. 
l2Q05l ). and four long exposures of 3060 s each were taken. 
The basic reduction tasks — trimming, bias and flat field 
correction — were done with the Gemini reduction package 
available for iraf. The extraction of the spectra was then 
performed using the task APALL under the National Optical 
Astronomy Observatories (NOAO) long-slit reduction pack- 
age. The spectra were extracted using the variance- weighted 
mode, with simultaneous extraction of the galaxy, sky and 
sigma spectra. The long GMOS slit (330 arcsec) allow us to 
measure the sky background in the same exposure as the 
galaxy spectrum, selected to be at the end of the slit. 

From the short exposure observations, we extracted the 
nuclear spectra with an aperture of diameter of 1.5 arcsec. 
From the long off-centre exposures, two apertures were ex- 
tracted in each exposure: one covering the radius from 30 
to 60 arcsec (hereafter, position 1), and another covering 
from 1 to 2 arcmin (hereafter, position 2). In Fig. [1] it is 
shown a Digital Sky Survey image of M32 with the GMOS 
slit indicated, positioned for the long exposures. The aper- 
tures extracted as positions 1 and 2 are indicated as filled 
regions in the slit, as well as the region used to measure the 
sky background. The signal from position 1 is on average 
four times stronger than the background level (from eight 
times the sky level on the side closest to the nucleus, to 
three times in the darker side). The signal from position 2 
is on average two times stronger than the background level 
(from three to 1.5 times). Position 2 was chos en to s ample 
the same galactic radius observed by Grillmair et al.l (|l996h 
in its resolved photometric study, whose approximate field 
of view is indicated in Fig. [T]as a square. 

The four spectra extracted at the central region have 
S/PsQ S/N ranging from 34 to 43 pixel"^ at 4750 A, and 
from 180 to 230 pixel" ^ at 7000 A. In the case of the long off- 
centre exposures, the spectra extracted for position 1 have 



Our observational data are described in Section 2, and 
we explain the analysis method in Section 3. Results are 
given in Section 4 and conclusions in Section 5. 



1 Computed as S/N= Sgalaxy/V^galaxy + Sgky, where Sgalaxy 

and 5sky are the galaxy and sky spectra, respectively. 
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S/N ranging from 50 to 80 and from 270 to 470 pixel" \ at 
4750 and 7000 A respectively. The spectra for position 2 have 
S/N from 27 to 46 and from 152 to 267 pixel" ^ at the same 
wavelengths. The S/N of the two off-centre apertures drop 
rapidly beyond 7000 A due to residuals of sky subtraction. 

We have not flux calibrated the spectra as no spec- 
trophotometric standard was observed in the nights of the 
observations. This will not hamper our analysis, as the full 
spectrum fitting technique adopted is not sensitive to the 
shape of the continuum. 



3 ANALYSIS 



Full spectrum fittin g techniq ues fe.g. jPanter et al. l2003l: 
Cid Fernandes et al.1 1 2005: M athis et al.ll2Q06l: lOcvirk et al.1 
2006al lbl: 'Walcher et al.i »2006i : iKoleva et al.ll2008l " have been 
used and validated in recent years as a powerful method for 
stellar population studies, and allows for fully exploiting the 
information present in our high-quality spectra of M32. 

We employed a mod i fied v ersion of the starlight code 
bv ICid Fernandes et al.l (|2005h in order to analyse the ob- 
served spectra. STARLIGHT combines spectra from a user- 
defined base of individual spectra in search for linear combi- 
nations that match an input observed spectrum. We adopt 
medium spectral-resolution simple SSP models spanning dif- 
ferent ages and metallicities as our spectral base (described 
further in the text), i.e., we describe the data in terms of a 
superposition of multiple bursts of star formation. 

To explore the parameter space defined by the linear 
coeficients associated with each population, the code uses a 
series of Markov Chains, each starting from different loca- 
tions. The likelyhood surface is explored with a simulated 
anneaUng scheme, consisting of a succession of steps where 
values are scaled by a gradually decreasing artificial 'tem- 
perature' (in analogy with Boltzmann factors in statistical 
mechanics), and parameter changes are performed with the 
Metropolis algorithm. In the computation of each spec- 
tral pixel is weighted by its correspondent error, but bad 
pixels (the two CCD gaps in CMOS), strong telluric absorp- 
tion features and residuals from sky subtraction are masked 
out. Pixels which deviate by more than three times the rms 
between the observed spectrum and the initial estimate for 
the model spectrum are also given zero weight. The inter- 
nal kinematics is determined simultaneously with the pop- 
ulation parameters. The fitted coefficients define a popula- 
tion vector (light fractions Xj at a reference wavelength of 
5000 A), which gives the combination of SSP models in the 
base that best fits (smallest x^) the observed spectrum. This 
same combination of techiques has proved useful in several 
studies of the sta r formation history of galaxies of different 
types (see, e.g., Asari et al.l l2007l and the starlight user 
manuao for examples) . 

For the present work, the public version of the code was 
modified to allow the use of continuum-normalized spectra 
(as opposed to fiux-calibrated spectra). A simple moving av- 
erage scheme of 200 A windows (excluding masked points) 
was used to remove the continuum shape of both model 



and observed spectra. This method fits all the pixels in- 
dependently of the continuum, being therefore insensitive 
to extinction or fiux calibration errors. The differences be- 
tween fiux-calibrated versus continuum-normaliz ed spectral 
fitting have not been extensively studied yet. I Wolf et al.l 
(|2007l ) compare four methods of deriving population pa- 
rameters from integrated stellar cluster spectra, and their 
results favour the use of continuum-normalized spectra, as 
this method performed better in simultaneously retrieving 
age and metallicity for the clusters. We also performed a pre- 
liminary test using both ve rsions of starlight on the public 
globular cluster spectra bv lSchiavon et al ] (|2005l l. We found 
that the two code versions yield results that agree typically 
within 0.2 dex in age, and 0.3 dex in metallicity. It is pos- 
sible, therefore, that the same analysis applied on a set of 
fiux-calibrated spectra would retrieve slightly different re- 
sults than the ones presented here. 

Spectral SSP models, adopted as base spectra, are un- 
dergoing constant revisions, following the improvement of 
their two main ingredients: stellar tracks and stellar spec- 
tral libraries. As our aim is to shed some light on the com- 
posite nature of M32 stellar population, without being bi- 
ased by a particular choice of models, we opted for per- 
forming the analysis using four different sets of models: 
BC03, PEGASE-HR, Chariot & Bruzual and Vazdekis et 
alB The two latter are preliminary versions of new mod- 
els being devel oped on the basis of the work by BC03 and 
IVazdekisl(|l999l l respectively, using the MILES (Medium res- 
olution INT Library of Empirical Spectra) empirica l library 
(jSanchez-Blazauez et al.]l2006al : ICenarro et ahlboOTl l. BC03, 
Chariot & Bruzual and PEGASE-HR models were com- 



puted adopting evolutionary tra cks bvlBertelli et al.| (| 19941 ). 
and Vazdekis adopted tracks bv ICirardi et al.l (|200ol ) . BC03 
and Chariot & Bruzual adopt initial mass function (IMF) 
from Chabrier (2003) and PEGASE-HR and Vazdekis et al. 
adopt the IMF from Salpeterj (|1955l ). For the fit procedure, 
we selected a range of ages between 100 Myr and 15 Gyr, 
spaced in log(age) = 0.25 dex. For the metallicites Z we 
adopted the full range provided by each set of models. The 
minimum wavelength of the fit is set by the observations 
(4700 A), and the maximum is set by the base models (see 
Table [U). 

The chosen models are all based on empirical stellar 
spectral libraries and as such are characterized by the abun- 
dance p attern typical of t he Milky Way. Models such as 
those bv lCoelho et al ] (j2002) can in principle be used with 
spectral fitting techniques to extract information of the a- 
enhancement in integrated light-studies. We do not consider 
such models in the present study, however, as we understand 
that an analysis of the a-enhancement would add an extra 
free parameter to an alread y delicate analy sis. This will not 
hamper our conclusions, as lworthevi (|2004l ) has shown that 
non-solar abundance ratio effects are mild in the case of M32, 
and do not change the population parameters significantly. 

In total 48 starlight fits were performed (three galac- 
tic positions x four exposures for each position x four sets of 
base models). Examples of the fits performed by starlight 
are illustrated in Figures O [3] and |4] for the models BC03, 
Pegase-HR and Vazdekis et al. respectively. A measure of 



http://www.starlight.ufsc.br/ 



^ http: / /www. iac.es/galeria/vazdekis/vazdekis_models_ssp_seds. html 
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Figure 1. Digitized Sky Survey image of the field containing M32. The GMOS sht positioned for the long-exposures is illustrated, with 
Position 1, Position 2 and sky background apertures indicated as filled reffl ons inside the sli t . The GMOS slit is in total 330 arcsec long. 
The square indicates the approximate position of the HST field studied by ICrillmair et al 1 (119961) . The two discontinuities indicated in 
the slit correspond to the CCD gaps. 



the quality of the fit is given by the average absolute devia- 
tion A: 



-y 



N 



(/model (A) — /obs(A)) I 
/obs(A) 



(1) 



where N is the number of pixels, /model is the fitted spec- 
trum and fobs is the observed spectrum. Typically the model 
spectra reproduces the observed within 2 per cent of residual 
flux (see Table [2]). 



4 RESULTS 

The analysis was performed with each observed spectrum 
individually, and the light-fraction vectors Xj obtained by 
STARLIGHT were averaged over the four exposures for each 
combination of base models and radius. The results of Xj 
are qualitatively shown in Fig. O as a function of log(age) 
and [Z/Zq], and in tabular form in the Appendix A. The 
projections of the averaged Xj in the log(age) and Z planes 
give the age distribution functions (ADFs) and metallicity 
distribution functions (MDFs), shown in Figs [6] and [3 re- 
spectively. 

It is worth remembering that the inversion from an ob- 
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Figure 6. Age distribution functions obtained for each observed 
aperture (in rows) and each set of models (in columns), according 
to the labels. 
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Figure 2. Observed spectra of M32 (black lines) and model spectra (red lines), for the three different observed apertures as indicated 
by the labels. Masked and rejected pixels are indicated as grey points. The average absolute deviation A of the fit (see text) is shown in 
each panel. Base models for these runs are BC03. 
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Figure 3. Same as in Figure [21 using PEGASE-HR as base models. 



served spectrum to a set of parameters is a complex and 
degenerate problem, and we should a void over detailed de - 
scriptions of the population mixture (|Ocvirk et al ] l2006al ). 
With that caveat in mind, we draw the following conclusions. 



ancient population were obtain ed from photome t ric stu dies, 
both for the central region (|Brown et al.1 |200Q| . l2008l ) and 
the out er reaches of the galaxy, ^ 3 arcmin away from the 
centre (jAlonso- Garcia et al.ll20oj V 



4.1 The distribution of ages 

The age distribution has two to three peaks, depending on 
the base model. Signal of an ancient population [log(age) 
^ 10] is retrieved in all positions. Previous evidence of the 



Intermediate- age populations [log(age) between 9.0 and 
9.9] are also always retrieved. Most of the previous spectral 
studies fin d a mean interm ediate age for the population of 
M32 (e.g. iRose et al.1 l2005l find a SSP age between 4 and 
7Gyr, depending on the radius), and the existence of such 
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Figure 4. Same as in Figure [2] using Vazdekis et al. as base models. 
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Figure 5. Representation of the light-fraction vectors Xj as a function of log(age) and [Z/Z©]. The results for the nucleus, positions 1 
and 2 are shown at the top, middle and bottom rows, respectively. Each column shows the results for a different set of base models, as 
indicated at the top of the figure. The light-fractions are given in tabular form in the appendix. 
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Figure 7. Metallicity distribution functions obtained for each 
observed aperture (in rows) and each set of models (in columns), 
according to the labels. The dashed line in the panels for Position 
2 correspond to the empirical metallicity distribution function by 
iGrillmair et alJ Ijldm) . 

popul ation also has been suppor ted by the photometric stud- 
ies bv lPavidge JensenI (|2007[ ) and references therein. 

In many cases, the inversion algorithm retrieves small 
amounts of a young population [log(age) < 9.0] ), typically 
stronger in the runs with PEGASE-HR models. iRose et aU 
§005) noted that a population in M32 cannot be any 
younger than 2 Gyr, or it would produce a strong signature 
in the Ca II index. As our spectrum does not extend blue- 
ward than 4800 A, the Ca II region is not considered in the 
fit. The determination of ages in integrated spectra may be a 
difficult problem, as extended horizontal branch morpholo- 
gies and/or blue stragglers, un accounted for in the SSP mod- 
els, can mimic youn^ a^es (Ide Freitas Pacheco Barbuvl 
I1995I : iMaraston ThomasI I2OOOI : iTrager et al.1 l2005l l. In- 
terpretation of integrated globular cluster spectra from 
ISchiavon et aD ([2005) with STECMAP (STElla r Content 
via Maximum A Posteriori; lOcvirk et al.[ l2006al ) supports 
the idea that blue HB stars in old metal-poor populations 
might ap pear as a spurio us young burst of star formation 
(|Koleva et al.l l2007l , I2OO8I , Ocvirk, private communication). 
Ho t HB stars in the ce ntral regions of M32 have been found 
by iBrown et al.l (|2000l l , and thus it is possible that a least 
a fraction of the young population retrieved here is a signal 
of these hot old stars. 

4.2 The distribution of abundances 

The MDFs functions are smoother than the age distribu- 
tions, with one or two peaks. In most of the cases, the peak 
of the metallicity distribution is centred at sub-solar values. 

Position 2, the farthest from the nucle us, samples th e 
same radius that was observed by Grillmai r et al.l (jl996[ ). 
The metallicity distribution obtained in that work is over- 
plotted to our distributions as a dashed line in Fig. [T] Our 
results are in agreement with the conclusion by Grillmair et 
al. that the metallicity distribution in M32 is broad. On the 



other hand, they obtain a distribution which is poor in stars 
with [Z/Zq] ^ —1. This is, in contrast, an important fraction 
in our mixtures regardless of the base models adopted. 

Several effects may be hampering the comparison be- 
tween the photometric a nd the spectroscopic me tallicity dis- 
tributions. The MDF bv lCrillmair et al.1 (|l996h is based on 
the colour distribution of red giant branch (RGB) stars, and 
as such it is affected by observational errors and age ef- 
fects (their distribution assumes a constant age of 10 Gyr). 
Furthermore, the set of underlying iso chrones of the spec- 
tral models used b y the present work (|Bertelli et al ] |l994l : 
ICirardi et al.ll200"ol l is different from the isochrones used fo r 
the star counts in Grillmair et al. (|Worthev et al.l [l994l V 
Moreover, the MDF by Grillmair et al. is in the form of 
number of stars per Z bin, i.e., N(Z)dZ, while ours is a 
light-fraction MDF, Ltot(Z)dZ. In a general case where Ltot 
is not constant with Z, then N{Z) and L{Z) are not the same, 
and one should not expect a perfect match. As pointed out 
by the anonymous referee, one might not expect an MDF 
derived from spectroscopic data to match an MDF derived 
from photometric data unless biases and observational errors 
are properly taken into account. 

In the mixtures retrieved by our analysis, the metal- 
poor component corresponds to a population younger than 
the metal-rich one (according to Fig. [5] and Tables Al— A9). 
It would be particularly puzzling if this trend is real, as it 
is in general expected that younger populations are formed 
from interstelar material already processed by older popu- 
lations, and thus should be more metal-rich. On the other 
hand, in some cases this younger metal-poor population ap- 
pears at the edges of the parameter space covered by the 
models, what might be an indication of an artefact effect. 

4.3 Mean population parameters 

The majority of the previous studies of M32 stellar pop- 
ulations publish the age and metallicities derived by the 
comparison to SSP models. Contrary to what is some- 
times assumed, the SSP- equivalent ages and metallicities 
are not the same as the luminosity- weighted age and 
metallicites (see e.g.lTrag;er et al. l l2000l : IWalch~ et al. I I2OO6I : 
iTrager Somervillell2009f L 

We computed the luminosity-weighted parameters of 
our population mixtures by adopting: 

(log (age)) = ^xj ■ log(age)^- (2) 

j 

and 

{[Z/ZQ]} = Y,^i-[Z/ZQh (3) 

j 

where the vector Xj gives the normalized light- fraction of 
the jth SSP component of the model fit — The 

SSP-equivalent parameters, on the other hand, are the age 
and metallicity of the jth SSP which best fits the observed 
spectrum (smallest A). 

We list both the luminosity-weighted and SSP- 
equivalent parameters in Table [21 derived independently 
from each of the four exposures. It can be seen that the 
difference from adopting distinct base models is larger than 
the variance of the starlight runs on the four exposures of 
each position. 
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Table 2. Light- weighted mean parameters and SSP- equivalent parameters obtained from the fits. Results are given for individual observed 
spectra. 



Light-weighted parameters Best SSP fit parameters 



Region 


Obs # 


Age[Gyr] 


[z/Zo] 


A[%] 


Age[Gyr] 


[z/Zoi 


A[%1 








BC03 










Nucleus 


1 


4.0 


0.12 


0.93 


3.2 


0.07 


0.94 




2 


4.5 


0.09 


0.98 


3.2 


0.07 


0.98 




3 


4.6 


0.07 


0.97 


3.2 


0.07 


0.97 




4 


4.1 


0.12 


1.00 


3.2 


0.07 


1.00 


Position 1 


1 


5.4 


-0.13 


1.14 


10.0 


-0.33 


1.18 




2 


3.7 


-0.38 


1.09 


15.0 


-0.33 


1.14 




3 


3.8 


-0.45 


1.14 


15.0 


-0.33 


1.20 




4 


3.6 


-0.49 


1.19 


15.0 


-0.33 


1.24 


Position 2 


1 


7.8 


-0.55 


2.11 


10.0 


-0.33 


2.13 




2 


4.6 


-0.65 


1.51 


5.5 


-0.33 


1.52 




3 


5.6 


-0.66 


1.54 


10.0 


-0.33 


1.57 




4 


4.6 


-0.63 


1.65 


10.0 


-0.33 


1.65 






Chariot & Bruzual 








Nucleus 


1 


6.1 


-0.05 


0.89 


5.5 


0.07 


0.95 




2 


6.4 


-0.03 


0.94 


5.5 


0.07 


0.99 




3 


6.8 


-0.02 


0.93 


5.5 


0.07 


0.96 




4 


6.3 


-0.04 


0.96 


5.5 


0.07 


1.00 


Position 1 


1 


6.8 


-0.21 


1.02 


10.0 


-0.33 


1.04 




2 


10.0 


-0.25 


1.05 


10.0 


-0.33 


1.09 




3 


11.7 


-0.26 


1.01 


10.0 


-0.33 


1.04 




4 


13.3 


-0.28 


1.03 


10.0 


-0.33 


1.05 


Position 2 


1 


11.5 


-0.49 


1.93 


10.0 


-0.33 


1.95 




2 


10.8 


-0.30 


1.33 


5.5 


-0.33 


1.35 




3 


12.6 


-0.33 


1.37 


10.0 


-0.33 


1.40 




4 


12.8 


-0.37 


1.45 


5.5 


-0.33 


1.49 








Pegase-HR 








Nucleus 


1 


1.7 


-0.03 


1.04 


3.0 


0.07 


1.04 




2 


1.5 


-0.04 


1.10 


3.0 


0.07 


1.11 




3 


2.0 


-0.03 


1.07 


3.0 


0.07 


1.07 




4 


1.9 


-0.03 


1.10 


3.0 


0.07 


1.10 


Position 1 


1 


3.6 


-0.37 


1.13 


6.0 


-0.33 


1.18 




2 


1.7 


-0.25 


1.11 


6.0 


-0.33 


1.13 




3 


1.9 


-0.25 


1.10 


6.0 


-0.33 


1.12 




4 


2.0 


-0.25 


1.13 


6.0 


-0.33 


1.15 


Position 2 


1 


8.8 


-0.78 


2.06 


3.0 


-0.33 


2.10 




2 


4.2 


-0.49 


1.23 


3.0 


-0.33 


1.27 




3 


4.7 


-0.52 


1.28 


3.0 


-0.33 


1.32 




4 


5.6 


-0.55 


1.31 


3.0 


-0.33 


1.35 








Vazdekis et al. 








Nucleus 


1 


2.3 


-0.02 


0.97 


3.2 


0.05 


0.97 




2 


2.4 


0.02 


1.01 


3.2 


0.05 


1.02 




3 


2.7 


-0.08 


1.01 


3.2 


0.05 


1.01 




4 


2.5 


-0.01 


1.04 


3.2 


0.05 


1.04 


Position 1 


1 


3.9 


-0.37 


1.07 


5.6 


-0.33 


1.11 




2 


4.5 


-0.30 


1.13 


3.2 


0.05 


1.17 




3 


5.0 


-0.30 


1.15 


3.2 


0.05 


1.18 




4 


5.3 


-0.38 


1.23 


5.6 


-0.33 


1.26 


Position 2 


1 


4.3 


-0.68 


1.95 


5.6 


-0.33 


1.99 




2 


5.2 


-0.53 


1.39 


3.2 


-0.33 


1.38 




3 


6.5 


-0.52 


1.41 


5.6 


-0.33 


1.44 




4 


8.0 


-0.58 


1.52 


10.0 


-0.63 


1.52 
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Table 3. Average luminosity- weighted parameters. 



Region 


Age (Gyr) 


[Z/Zo] 


Nucleus 
Position 1 
Position 2 


BC03 
4.3 ± 0.3 
4.1 ± 0.9 
5.7 ± 1.5 


0.1 ± 0.1 
-0.4 ± 0.2 
-0.6 ± 0.1 


Chariot & Bruzual 
Nucleus 6.4 ± 0.3 0.0 ± 0.1 
Position 1 10.4 ± 2.8 -0.2 ± 0.1 
Position 2 11.9 ± 0.9 -0.4 ± 0.1 


Nucleus 
Position 1 
Position 2 


Pegase-HR 
1.8 ± 0.2 
2.3 ± 0.9 
5.8 ± 2.1 


0.0 ± 0.1 
-0.3 ± 0.1 
-0.6 ± 0.1 


Nucleus 
Position 1 
Position 2 


Vazdekis et al. 
2.5 ± 0.2 
4.7 ± 0.6 
6.0 ± 1.6 


0.0 ± 0.1 
-0.3 ± 0.1 
-0.6 ± 0.1 



4.4 Evidence for gradients ? 

Studies of radial colour and line strength gradients have 
led to disc ordant results in recent years (see discussion in 
iRose et al] [2005. and references therein). In the present 
work, we opted for maximising th e S/N (a criti cal issue in 
spectral fitting technique, see e.g. iKoleva et all 2007^. and 
thus we have not split our spectra in many apertures to 
perform a detailed gradient study. In general terms, it is 
possible to identify in our analysis a shift in the distribution 
of populations, in the sense that the fraction coming from 
an old population is higher in the outer regions than in the 
nucleus. In order to statistically address this issue, in Table 
[3] we show the light- weigthed parameters from Table [2] av- 
eraged on the four exposures, with their corresponding 1 a 
uncertainties. A trend of higher ages and lower metallicities 
as one moves from the nucleus to position 2 can be seen in all 
cases, even though the gradient slopes are model-dependent. 

These fi ndings supports one of the scenarios discussed 
bv iRose et a l. (2005), in which there is a shifting balance 
between the populations with radii. A recent burst of star- 
formation on top of an ancient population would appear as 
a higher concentration of a younger and more metal-rich 
population in the nucleus. This scenar io would be a n atu- 
ral result of the model proposed by Bekki et al.l (|200lh , in 
which centralised star formation in M32 is triggered by the 
tidal field of M31. On the other hand, HST ne ar- ultraviolet 
(UV) imaging of M32 by iBrown et aP (jlQQSl ) reveals that 
the resolved UV-bright stars are more centrally concentrated 
than the underlying unresolved UV sources. Thus, it is also 
possible that the trend we see is a changing influence of a 
relatively small hot population prevalent in the nucleus. As 
the models do not yet consider hot stars in old populations, 
it is not possible with the present analysis to distinguish if it 
is the hot-old stars or the hot-young stars (or a combination 
of both) which is driving this balance shift 

^ Note added in proof: An alternative scenario to explain our 



4.5 Dependence of the results with wavelength 
range 

For the main fitting runs we opted for maximizing the use 
of information and thus to perform the fitting on the to- 
tal wavelength range allowed by each set of spectral mod- 
els. We additionally performed tests fitting only the region 
which is common to all base models (4700-6800 A) and 
a limited region encompassing the most commonly used 
Lick/IDS (Lick Observatory Image Dissector Scanner) in- 
dices (4828-5363 A, from H/3 to Fe5335). 

The results fitting only the region 4700—6800 A are not 
strongly different from those obtained when fitting the whole 
range, except that lower weight is given to the old popula- 
tions, and thus mean ages which are roughly 1 Gyr younger 
are derived. In the case of fitting the 4828— 5363 A re- 
gion, the ages retrieved are slightly older in the nucleus 
but younger in the outer regions, and thus the gradient 
is flatter. The MDFs do not change significantly, except 
that the mean metallicity is slightly higher when only the 
4828-5363 A range is fitted. 

The dependence of the results on wavelength range is 
an expected behavior for galaxies with star formation his- 
tories that are not well represe nted by a singl e burst (e.g. 
ISanche z-BlazQ uez et all r2QQ6bl : ISchiavonI l2007f ) . Therefore, 
we favor the use of wide wavelength ranges in analyses that 
aim at studying composite stellar populations. 



5 SUMMARY AND CONCLUSIONS 

We obtained high S/N integrated spectra of the galaxy M32 
with the GMOS at Gemini-North, out to 2 arcmin (3.6reff) 
from the centre of the galaxy. We performed an analysis of 
the composite stellar population in M32 using starlight 
by Cid Fernandes et al. (2005), a technique of full spectrum 
fitting, and adopting different set s of medium spectral reso - 
lution stellar populatio n models: IBruzual CharlotI (|2003l ), 
iLe Borgne et al.l (|2004h . Chariot & Bruzual (in prep.) and 
Vazdekis et al. (in prep.). Three different regions of the 
galaxy were analysed separately, the nuclear region, a re- 
gion from 30 to 60arcsec, and a region from 1 to 2 arcmin. 

We can find mixtures of populations that reproduce our 
observed spectra within 2 per cent of residual fiux. We find 
that in the nuclear region, 23—65 per cent (depending on the 
model used) of the light at 5000 A comes from an ancient 
population [older than log(age) = 10], and intermediate- 
age populations [log(age) between 9 and 9.9] contribute 
to 25— 48 per cent of the light. Beyond 1 reff (reff ~ 33"), 
the light coming from the ancient population increases 
(27—90 per cent) and the contribution from intermediate- 
age populations decreases (4—41 per cent). There is a signal 
of a young-population [log(age) younger than 9.0] that con- 
tributes up to 30 per cent. As no evidence of such young 

findings is the one proposed recently by Kormendy et al. (2009). 
In this formation scenario, the young light in the centre would 
be the result of a starburst in the latest dissipative merger that 
made the galaxy. There is a gradual change as radius increases, 
towards larger contributions from older stars that formed before 
the most recent merger (s). In this scenario, the gradient would 
be therefore intrinsic to the formation process of M32, and not 
triggered by M31. 
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population has been found by previous photometric and 
spectroscopic studies, it is hkely that at least part of this 
populati ons is a spect ral signal of the hot HB stars iden- 
tified by iBrown et al ] liooo). These findings suggest that 
the lack of spectral population models accouting for these 
hot old stars may set a constraint in our ability to identify 
the detailed age distribution of a galaxy from its integrated 
spectrum. 

The metallicity distribution is peaked at subsolar val- 
ues, and we retrieve a non-negligible contribution from a 
metal-poor population with [Z/Z©] ^ —1.0. Such metal- 
poor p opulation has no t been identified in the CMD analysis 
such as lGrillmair et aP (|l996l ). As our position 2 samples the 
same radius as the analysis by Grillmair et al., this disagree- 
ment between integrated light studies and resolved stellar 
studies has yet to be understood. 

As the main result of the present analysis, we propose 
that an ancient and an intermediate- age population co-exist 
in M3 2, in agreement with the photometric evidence f ound 
by e.g. iBrown et aP (|200Q|1: I Alonso- Garcia et all (|2004l ) (for 
an old population) and Davidge & Jensen (|20Q7l ) (for an in- 
termediate age population). Moreover, we propose that the 
balance between these two populations change from the nu- 
cleus to the halo (outside Ireff) in the sense that the con- 
tribution from the intermediate population is larger at the 
nuclear region. 

More detailed descriptions of the age and metallicity 
distributions in M32 may have to wait for future stellar 
population models, as the balance of populations that fit 
our high-quality integrated spectra is model-dependent. 
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APPENDIX A: RESULTS 

Population mixtures obtained for each set of SSP models, 
averaged on the four exposures (Tables Al— A12). 
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Table Al. Light percentage fraction as a function of age and Z 
for the nuclear region, obtained with BC03 models. The values 
for [Z/Zq] are given in the column headers, log(age) is given in 
the first column. Adopted: Zq = 0.017. 



[Z/H] 
log(age) 


-2.23 


-1.63 


-0.63 


-0.33 


0.07 


0.47 


8.00 
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36 





9.25 

















5 


9.50 




















9.75 














3 





10.00 














38 
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9 
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Table A2. Same as in lAll but for position 1 and BC03 models. 



[Z/H] 
log(age) 


-2.23 


-1.63 


-0.63 


-0.33 


0.07 


0.47 


8.00 
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Table A3. Same as in lAll but for position 2 and BC03 models. 



[Z/H] 
log(age) 


-2.23 


-1.63 


-0.63 


-0.33 


0.07 


0.47 


8.00 
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19 
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6 
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Table A4. Same as in lAll but for the nuclear region and Chariot 
&; Bruzual models. 



[Z/H] 
log(age) 


-2.23 


-1.63 


-0.63 


-0.33 


0.07 


0.47 


8.00 
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Table A5. Same as in lAll but for position 1 and Chariot & 
Bruzual models. 



[Z/H] 
log(age) 


-2.23 


-1.63 


-0.63 


-0.33 


0.07 
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8.00 
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1 
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Table A6. Same as in lAll but for position 2 and Chariot &: 
Bruzual models. 



[Z/H] 
log(age) 


-2.23 


-1.63 


-0.63 


-0.33 


0.07 
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Table A7. Same as in lAll but for the nuclear region and 
PEGASE-HR models. 



Table A 10. Same as in lAll but for the nuclear region and 
Vazdekis et al. models. 



[Z/H] 
log(age) 


-1.63 


-0.63 


-0.33 


0.07 


0.47 


[Z/H] 
log(age) 
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-1.23 


-0.63 


-0.33 


0.05 


0.25 


8.00 





29 











8.00 








3 





13 





8.25 

















8.25 


7 

















8.50 

















8.50 




















8.75 

















8.75 




















9.00 

















9.00 




















9.25 

















9.25 




















9.50 














25 


9.50 

















43 


9.75 














23 


9.75 




















10.00 








23 








10.00 




















10.18 

















10.18 











4 


29 






Table A8. Same as in|AT] but for position 1 and PEGASE-HR 
models. 



Table All. Same as in lAll but for position 1 and Vazdekis et 
al. models. 



[Z/H] 
log(age) 


-1.63 


-0.63 


-0.33 


0.07 


0.47 


[Z/H] 
log(age) 


-1.63 


-1.23 


-0.63 


-0.33 


0.05 


0.25 


8.00 





28 


3 








8.00 





6 














8.25 

















8.25 


6 








5 








8.50 

















8.50 




















8.75 

















8.75 




















9.00 

















9.00 




















9.25 

















9.25 











1 








9.50 

















9.50 

















4 


9.75 








17 


6 


22 


9.75 














36 





10.00 


6 





21 








10.00 




















10.18 

















10.18 





15 








27 






Table A9. Same as in|AT] but for position 2 and PEGASE-HR 
models. 



Table A12. Same as in lAll but for position 2 and Vazdekis et 
al. models. 



[Z/H] 
log(age) 


-1.63 


-0.63 


-0.33 


0.07 


0.47 


[Z/H] 
log(age) 


-1.63 


-1.23 


-0.63 


-0.33 


0.05 


0.25 


8.00 





11 











8.00 





5 














8.25 

















8.25 


14 








1 








8.50 

















8.50 




















8.75 

















8.75 




















9.00 

















9.00 




















9.25 

















9.25 




















9.50 

















9.50 




















9.75 














17 


9.75 











4 


5 





10.00 


25 


24 


19 





4 


10.00 




















10.18 

















10.18 





25 








47 






